V in other coordinates 1 


Chapter 15 V in other coordinates 


On a number of occasions we have noticed that del is geometrically determined — it does 
not depend on a choice of coordinates for IR”. This was shown to be true for Vf, the gradient 
of a function from R” to R (Section 2H). It was also verified for V e F’, the divergence of a 
function from R” to R” (Section 14B). And in the case n = 3, we saw in Section 13G that it 
is true for V x F, the curl of a function from R° to R?. 

These three instances beg the question of how we might express V in other coordinate 
systems for R". A recent example of this is found in Section 13G, where a formula is given for 
V x F in terms of an arbitrary right-handed orthonormal frame for R°. We shall accomplish 
much more in this chapter. 

A very interesting book about V, by Harry Moritz Schey, has the interesting title Div, 
Grad, Curl, And All That. 





















































A. Biorthogonal systems 


We begin with some elementary linear algebra. Consider an arbitrary frame {y1, Y2,-.-, Yn} 
for R". We know of course that the Gram matrix is of great interest: 














G = (gi © Yj). 
This is a symmetric positive definite matrix, and we shall denote its entries as 
A oe at dea oe 


Of course, G is the identity matrix <=> we have an orthonormal frame. 
We also form the matrix ® whose columns are the vectors y;. Symbolically we write 


® = (Y1 Y2--- Yn): 


We know that ® is an orthogonal matrix <=> we have an orthonormal frame (Problem 4-20). 
Now denote by WV the transpose of the inverse of ®, and express this new matrix in terms 
of its columns as 


W = (dy We... Un). 
We then have the matrix product 
vy 
Wd =] : | (y1...@n) 
yt 
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But since W equals the inverse of ®, we conclude that 


Wi © Yj = Oyj. 


DEFINITION. Two frames {y1,...,¢%n} and {v1,..., an} are called a biorthogonal system 
if 


Wi @ pj = O53. 
Clearly, {¥1,...,n} and {Y1,...,¢n,} form a biorthogonal system <=> the frame {y,..., Yn} 
is an orthonormal one. So the present definition should be viewed as a generalization of the 
concept of orthonormal basis. 


Also it is clear that the relation between the y;’s and the y;’s given here is completely 
symmetric. 


PROBLEM 15-1. Let a frame for R? be {7, at + bj}, where of course b 4 0. Compute 
the corresponding {w, W2} which produces a biorthogonal system. Sketch all four vectors 
on a copy of R?. 


PROBLEM 15-2. Let {y1, 2,93} be a frame for R*®. Prove that the biorthogonal 
frame is given by 
1 * Pa 
by = PEM Pa 


etc. 
ly, p2; ~3] 


PROBLEM 15-3. Suppose {y1,...,Qn} is an orthogonal frame. Show that the 
corresponding vectors w1,..., Wn are given as 


Pi 
i= oo. 
Ilvall? 


PROBLEM 15-4. Given a frame {yj,...,¢,} denote 
J =det ®, 


Prove that det G = J”. 
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Finally, a little more notation. It has become customary to denote the inverse of G' = (g;;) 
by : 
Gt = (g"") 
(called “raising the indices”). That is, 


So on = > ging = ij. 
k=1 k=1 


Inasmuch as 


G='d, 
we conclude that 
U = (6')! 
= 6G"!. 


In terms of entries of these matrices this equation means that 


n 


(Wij = S-(®)ing”. 


k=1 
Thus we find that the columns satisfy 


n 


v= S org”. 


k=1 


Rewriting this equation gives 
Wi= ye 903. 
j=l 


The corresponding inverse equation is of course 


2 » Gj Vj: 
j=l 





PROBLEM 15-5. Prove that for all v € R” 











odes 
i=1 
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The result of this problem is that when you wish to express v as a linear combination of the 
basis vectors y;, the corresponding coefficients are given directly in terms of inner products 
with the y,;’s. You might think of the formal expression 


n 


Sidi e ())¢i 


i=1 





as representing the identity linear function on R”. In other words, 


PROBLEM 15-6. Prove that 


T= Ss pi}. 
i=1 


B. The gradient 


We continue with an arbitrary biorthogonal system on R”, and maintain the notation of 
Section A. 

Suppose that IR” _, R is differentiable at a point « € R”. Our task is to express the 
vector V f(x) € R” as a linear combination of the vectors y; in the given frame. This is quite 
an easy task. In fact, Problem 15-5 gives immediately 


























Via)= De Vi (2) © divi. 


The inner products in this formula are of course directional derivatives of f in the directions 
w;. Thus, in the notation of Section 2C 


Vilg)= >» PF le; Vi) Gi. 


PROBLEM 15-7. Show that everything in this formula can be expressed entirely in 
terms of the frame {y1,...,n} by a double sum 


Vi (2) = S° og DF (a; 9;)e- 


t,j=1 
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C. The divergence 





We want to discuss a vector field f defined on an open subset of IR”. We can thus regard 
f as a function from R” to R”, and as such it has a derivative. At a point x in its domain, the 
derivative Df (x) is a linear transformation of R” to R", represented in terms of the standard 
coordinate basis €1,...,é€,, by the n x n Jacobian matrix 


Of; 
(se). 


It is the trace of this matrix which is V e f, the divergence of f. This observation is the key 
to our representation of V e f, and we need a simple fact from linear algebra: 
































THEOREM. Suppose A is a real n x n matrix, and regard the y;’s and w;’s as column 
vectors. Then 


traceA = > Ad; © Yj. 


i=1 
PROOF. We use the matrices ® and W from Section A, so that 
Ay, ep; = ji entry of the matrix ® “AW. 
Thus 
3 Ay; ey; = trace(®*AV) 
~ = trace(AW®‘) (Problem 4-2) 
= trace(A), 


since W and © are inverses. 
QED 
As a result we now have 





Ve f(z) => 2h 


= OG, 
i=1 


=S 0 Df (a)bie ¢ 


i=1 


= > Df (x; Yi) © Yi. 


i=1 
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D. The curl 


Before going into the representation of curl, we summarize what we have obtained to this 
point. In terms of a fixed biorthogonal system we have found 


Vig)=Vi@)= > Pfla: Vi) Vi 


and 


divF(x) = Ve F(z) = es DF(a; Ui) © Yj. 


These two formulas can be incorporated as a single expression 


y= S° vi times D( ; uv) 


i=1 


where V and the right side are both regarded as operators. They operate on scalar-valued 
functions to produce the gradient and on vector fields to produce the divergence. In the 
gradient situation “times” is scalar multiplication, whereas in the divergence situation it is 
dot product. 

A good way to remember this formula is to replace the directional derivative 


Df (x; Vi) 


by the symbol 
Of 
Ov; 





Then we have 








” O 
v= 2 times ri 








This leads us to an educated guess for curl, where we just use the same formula. Thus we 
suppose that {y1, v2, 3} and {v1, v2, v3} form a biorthogonal system for R®. Then we have 





THEOREM. For a vector field F on R®, 





curlF (x) = V x F(x) = 50 gi x DF(x;y). 


i=1 


V in other coordinates 7 


PROOF. Every vector field F' can be expressed in the given frame in the form 


F(a) = > file). 


Thus it suffices to prove the formula for the special case of vector fields of the form f;(x)p;, 
for 7 = 1, 2, 3. Thus we assume 


F(x) = f(x)y;. 
We have from the product rule of Problem 13-14, 
V x F(x) =Vf (x) x vy; (since y; is constant) 
and our formula for gradient gives 


3 
V x F(2) = S° Df (a; videi ec) 


i=1 


3 
a Ye x Df (x; Wi) y; 


i=1 


= So ¢i x D( fp;) (a; Wi) 


i=l 
3 

a Ye x DF (x; Yi). 
i=1 


QED 


PROBLEM 15-8. Are you surprised that the formula for curl does not require the 
frame {(~1, Y2, 93} to be right-handed? Explain what happens if v3 is replaced with —y3. 


PROBLEM 15-9. Prove that 


Y1 X Yo = JY, 
po X v3 = JY, 
3 X Yi = Jy. 
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PROBLEM 15-10. Show how the formulas of the preceding problem behave when 
each y; is replaced with the vector a;y;, where the a;’s are nonzero scalars. 


E. Curvilinear coordinates 





All that we have done up to now is represent V in terms of a fixed frame for R”. In 
practice, however, the coordinates themselves are changed in a perhaps nonlinear fashion. We 
still need convenient expressions for V in terms of the new coordinates. 

We shall call our coordinate system T. What we mean is that T is a C! function from R” 
to R”, and has a C' inverse as well. (Usually the domains of T and T~! will not be all of R”.) 
The notation we shall use is 











z= T(u), 
representing the formula which determines x = (x1,...,2,) in terms of the curvilinear coor- 
dinates u = (u1,...,Un). The inverse function 

“e=T (x) 


gives the curvilinear coordinates of the given point x. 
The Jacobian matrix of T’ is of course 


Ox. eLae 
Ou, i Oun 
/ = ; : 
Ttuy=] : 
Orn Orn 
Ou, aoe Oun 


This matrix must be invertible and we denote its determinant by 


J(u) = det T’(u). 





The columns of 7’(u) represent linearly independent vectors of IR”, and we think of them as 
“attached” to the point x = T(u). That is, they are tangent vectors to the manifold R” at 
the point x. These vectors play the role of a moving frame for R”: 








OT 
C= Fan Lead 





There are other natural tangent vectors at x, namely the gradients of the cordinates u; viewed 
as functions of x. We denote them as 


w(x) = Vuy, Le 1. 
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Since u; is the j coordinate of T~!(x), the gradient Vu; is the j™ row of the Jacobian matrix 
for T~!. Therefore, since the chain rule implies 


(T)'(L(u))T'(u) = L, 


we have the matrix equation 


v1 
We 
| (Yr G2--- Pn) =I 
Un 
That is, 
Wie pj = ig. 


In other words, the frames at x given by 


{ei(T(@)),---,Pn(P(a@))F and {vi(x),.-.,Un(a)} 


are biorthogonal. Here’s a representative sketch in the x;—®2 plane: 


10 








Our point of view is going to be that when we work on 


U1,-- 
and the “natural” derivatives 0/Ou;. 


Of course, each 7; is a vector field on 
also a vector field on 
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= constant 
U,= 


= constant 
Uy= 


The Y ‘sare 
tangent to 

the coordinate 
curves 


The wy, ‘sare 
orthogonal to 
the coordinate 
curves 














R” and use the coordinates 


.,Un, We are interested in calculating everything in terms of the “natural” vectors y;(w) 














R” — in fact, a gradient field. And each y; 0 T is 


R”, but not necessarily a gradient field. 
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We still maintain the notations introduced in Section A. Thus 
gig(u) = vilu) © pj(u), 
and the matrix inverse to (g;;) has entries 
g"'(u). 
Furthermore, it follows that, with « = T(u), 


n 


bi(z) = >> 9 (uy; (u) 


and 


yi(u) = > Gi(u)e;(2). 


F. The gradient 
The formula for Vf goes over with no change from Section B. We read it off from Problem 
15-7: with z = T(wu), 
Vi (x) = d¢ 9 (u) DF (a; 9;(u)) piu). 
ij=l 


We notice that the chain rule gives 





Df (2; p;(u)) = VF(T(u)) © oT 
O 
= Te (T(u)). 


NOTATION. In this context we denote the pull-back of f by the coordinate transformation 
x = T(u) to be the function 


Thus we have 





vile)= >> (uw) sre) 


i,j=l 
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This is exactly the sort of formula we want. Everything on the right side of the equation is in 
terms of the coordinates we are working with. 


G. Spherical coordinates 





We pause for a significant example. We use our standard spherical coordinates for R?: 
x=rsinycos#, 
y=rsingsin#g, 


Z=T cosy. 
We then have the frame 


yi = (sin y cos 8, sin ysin 6, cos y), 
2 = r(cos ycos 8, cos y sin 6, — sin y), 
3 = rsin y(— sin 8, cos 8,0). 
Thus gi; = 0 if i Aj, and 
gu = 1, 
92 =T, 
933 = T’ sin” yp. 


Furthermore, g’? = 0 if i 4 j, and 


g? =1, 
1 
22 
g = pe 
33 1 
r2 sin? py 
Thus 
Of* Lay 1 «oF 


| | 

Or Ph" pp dy ae. sin? » 00 ve. 

Frequently this formula is displayed in terms of the orthonormal frame associated with {y1, Y2, y3}. 
Namely, define 





Vi(x,Y, 2) = 


~ _ Pi 
P= : 
Ilyall 





Then it becomes aft 
Vf= a git 
r 


1of* . 1 Of*. 
a4 
r Oy 





“ rsingy 06 a 
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In fact, a nice alternative for this notation is to say that in general 





ar 

“a Ou; 
Ui = Tor 
ll Sas 











so the formula becomes 





oT. 
r- 


TOF" A Of" 5 
Or r dp” 


Mi rsiny 00 — 


PROBLEM 15-11. Show that J =r’siny. 


PROBLEM 15-12. Carry out the same analysis for the case of cylindrical coordinates 

















for R?: 
tT =7 C0s0, 
y=rTemn7g, 
Z= 2. 
Conclude that ve Lae or 
vpn es Pp, OF. 


Or r 00...” @z 


H. The divergence 


A new situation arises in the curvilinear representation of divergence and curl. Namely we 
assume from the start that the vector field in question is represented in terms of curvilinear 
coordinates and the frame {y ,...,%,} associated with them. Thus we treat the vector field 
F(a) by first expressing it in the form 


fa) = d Filu)yi(u), where x = T(u). 


We then want to express the scalar Ve F' in terms of derivatives 0/Ou; involving the coefficients 
F;(u). Since the vector y;(u) comes from T by differentiation and is liable to be subject to 
differentiation again, we need to assume T is of class C?. We shall make this assumption for 
the remainder of the chapter. 

It turns out that there is a huge simplification available if we first establish the remarkable 
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LEMMA. Ve (J-!y;) = 0. 


PROOF. The vector field we are dealing with in this context is 
G(r) = J~*(u)y;(u), where x =T(u). 


We give a proof of this lemma based on the divergence theorem. To set this up, suppose h is 
any C' real-valued function on R” which is zero outside a small ball. We shall prove that 





[ rvecur=0 (*) 


Since V eG is a continuous function, this will prove the result. For if V e G(xo) > 0 for some 
Xo, then V e G(x) > 0 for all x in some neighborhood of xp. We can then choose a suitable h 
to make the integrand in (*) always > 0 and positive in a neighborhood of zo, so the integral 
in («) will be positive, a contradiction. Likewise if V e G(ao) < 0. 

The function A in this calculation is frequently called a test function. It is unimportant in 
itself, but is used to “test” the crucial function V eG. 

We now turn to the proof of (*). First, 


[ vetnayar =0 


follows from the divergence theorem. Actually, just the fundamental theorem of calculus is 
needed, as a typical term in the Cartesian formula for divergence is 


O 
—(hG eé;)dzx, 
aes Bx, | i) 


and performing the x; integration first gives the result, since h is zero outside a small ball. 
Since the product rule gives 


Ve(hG) =VheGthVeG, 


the proof of (*) now reduces to proving that 


VheGdz = 0. 


R” 
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In this integral we change variables with x = T(u), and we obtain 


| VheGdxr = | Vh(T(u)) e G(T(u))| det T’(u)|du 


R” 


VA(T(u)) ¢ J*(u)yi(u)|J(u)|du 


R” 


= es Vh(T(u)) © - 


a 


du 





chain rule O 


= (h(T(u))du 





FTC 
= 0. 


QED 
Now we compute the divergence of a typical summand in F'(x): the product rule gives 


V0 (F(u)pi(u)) = Ve JMR) Fr eilu) 


= V(I(u)F(u)) © 








Tray Pi): 
the other term being zero thanks to the lemma. And now we invoke the chain rule to get 
1 1 oT 
ae Vd (OE i(u) = ~V( SF; 
Fay Hw) Flu) « valu) = FVUIP) 0 5 
10 


Thus we have our final formula: in case 


Fa) = Awe(u), 2= Tu), 


then 





1 Qa J(u)F) 
VeF= : 
J(u) Ou; 
Notice how closely this resembles the Cartesian case! 
For our spherical coordinate example, write 














A 


EGY.) = Fi(r, p, 8)? + Foe + F308, 
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so that the corresponding functions in our formula are actually 


Thus 


PRO 


F,, Fo/r, and F3/rsiny, respectively. 








ae koe Cig. 0 
Vere ose ra sin pF) + ae oe + 59 (Fs) 
Par) + ingap eel) * rang a0 














BLEM 15-13. For the case of cylindrical coordinates for R® let 


F(a, y,z) = Fyf + Fo + Fé, 


and prove that 


PRO 


_ 10(rF\) 1 OF, OF: 
Me ae Orr OO Oz 














BLEM 15-14. Suppose that {y1,..., yn} is a frame field for R", so that it may 





vary from point to point. Suppose that the corresponding frame field {2 ,...,w,} yields 
a biorthogonal system. Prove that the basic formula of Section C is still valid: 


- OF 
VeF= pie =. 
DP og 


I. The curl 


We continue the notation of the preceding section, so that we are concerned with a vector 


field on | 





R®, but instead of expanding F in terms of the frame {¥1, ¢2, y3}, it seems better to 


work with F(x) itself. We then follow the notation of Section F and define the pull-back 


F*(u) = F(T(u)). 


We shall employ the formula for curl as given in Section D, but with the y;’s and y,’s inter- 
changed. So we find 


3 
V x F(x) =) >i x DF (2; )). 
4=1 
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Since y; = OT /Ou;, the chain rule gives with x = T(w), 





ee 
V x F(a = ows x =—(u). 


Next we insert the algebra of Problem 15-8: 
pi = JI "v2 x v3 ete, 
obtaining 


* 


Ou, 


The formula for vector triple product of Problem 7—4 gives 


Vx F=J7"(we x 3) x 





* * 


O 
F= ay, eee 
JV Xx p2e on p3 — y3e cn 





Yot::: 


The product rule gives 








O O 
IV x F=5-(peF)ps— 5 0 F's 
g) Ops 
ye F* F* bee: 
Du, (¢3 @ F™) po + co as 


The other eight terms are obtained from these by cycling through the indices 1-~2—3-—1. 
Look at the terms in which F* appears undifferentiated: the coefficient of y3 in these terms is 


But 
Ovi _ O OT 0 OT _ Ope 


Ou2 - Ouse Out ~ Oty dua Ou, : 
so that the coefficient of y3 just obtained is zero. Notice here that we have the assumption 
that T is of class C? so that the mixed partial derivatives are equal. Likewise for y, and yo. 
Thus we have 





O O 
JIVxF= pun ee a ee ee 
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(four more terms). This expression is precisely the formal determinant 
Y1 2 3 
O O 
det ? 
Ou, Ouse Ous 
piel” poek™ y3ek* 
This gives our final formula for curl: 
Y1 2 3 
A O O O 
V x F = —det 
J Out Juz Juz 
pie poek™ y3ef* 











Of course, the pattern is just that of the original definition of curl in terms of Cartesian 


A 


coordinates. 
Now let’s examine the spherical coordinate case: 


F(a, ; z) = Fi(r, yp, 0)F =F Pye + F306. 


Then 
pie F* =feh*=F, 
poe F* =rgpek* =rFy, 
3 @ F* =rsingé e F* =rsinyFs, 
so that 
r rp rsingé 


0 O 0 
06 





FP, rfy rsingFs 
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PROBLEM 15-15. For the case of cylindrical coordinates in the notation of Problem 


A 


15-12, show that 
r rd ¢Z 
Yxp=oat|2. 2. 2. 
r r OO Oz 
F TP» F, 


In particular, if the vector field is planar, 
P= Pylr, O)r + Fair, 0)0, 


conclude that 





_f era, 4 10F; \ > 
vxr=(F a a Nk 


It is especially interesting to see the formula for the operator V in terms of these curvilinear 


coordinates. The directional derivative in the direction y;(u), as we have noticed, is just 0/Ou;. 


Thus we may rewrite the general formula from Section D in the form 


Y= Soi times =. 
i=1 


or if we replace w by its definition, 











V= » Vu; times — 


i=1 








As always, times means 
e scalar multiplication if we are calculating gradient, 
e dot product if we are calculating divergence, 


e cross product if we are calculating curl. 
This version of the expression for V is rather easy to remember. If we think of V loosely 


as 0/Ox, then the formula is similar to 
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or, even more briefly, 


QO Ou ‘g 

= == limes —. 

Ox Ox Ou 
J. Orthogonal curvilinear coordinates 


Many coordinate systems that arise in specific applications have the feature that the curves 
in R” determined by them are orthogonal. This means in our notation that for i 4 7 


OT OT 0 
Pie Pj = ae = 
Ou; Ou; 
In other words, g;; = 0 if i A 7. Then the numbers of importance are 


Gn), Las, 


2 0 if i Xj, 
gtuy= 4 , © gee 
Gulu) 1 L= j. 








and 


Moreover, 
J? = 911922 - ++ Gnn- 
Notably, spherical coordinates and also cylindrical coordinates both fit this situation. 
Our formulas for V all simplify under the assumption of orthogonality. Thus the gradient 
is given by 


n 


Vi@= 2 Pa 


The divergence is still the same, no simplification occuring at all. 
The formula for curl is much simpler. For in the notation 


1 oOf* 





we now have simply 
pie FY = gif. 





Thus 
Y1 ~2 ~3 
Vx F= = det e d a 
911922933 Ou, Ou2 Ouz 


gl goek2 933h3 
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PROBLEM 15-16. If you use instead the associated orthonormal frame so that 


show that 





Ju: JGa22 /933P3 


1 O O 0 
V x F = —— det 
911922933 Ouy Ou2 Ou 








Jul VSg22F2 ./933F3 


kK. The Laplacian 


Without a doubt the most important second order partial differential operator on R” is 
the Laplace operator, or the Laplacian. It operates on a real-valued function by first forming 
its gradient and then forming the divergence of the resulting vector field. Thus the notation 


is V?: 





Vf =Ve(Vf). 


Another common notation is A = V?. 


Thus we have in Cartesian coordinates 


Of course, we well recognize that V? has an intrinsic geometric meaning apart from any 
coordinate system, since the same is true of Vf and of the divergence of a vector field. 


Before giving the general formula in curvilinear coordinates we mention the linear case in 
which x = T(u) is given by matrix multiplication 


x= Au, 
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where A = (a;;) is a real nonsingular n x n matrix. Then our notation gives the following: 


aij 
Ox a2; 
y,=~—=| ” | =j" columnof A; 
Ou; 
anj 
Wi = Vuj= V(Ao tx) en entry 
il 
i2 
= = transpose of i‘ row of A7?; 
qit 
J = det A; 


nm 

Gij = y Akidk;; 
k=1 
nm 

g! _ . ata, 
k=1 


PROBLEM 15-17. Prove that 





7 Ff) 
V?f (x) 2 Babes 


PROBLEM 15-18. Prove that V? is invariant with respect to this change of variables, 
that is, 

n ay” 

2 = 


i=l 7 





if and only if A € O(n). 
Now we turn to the general case of curvilinear coordinates. We have from Section F 
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so the gradient vector field V f is represented as in Section H with the corresponding coefficients 
denoted there as F;(u) given by 


Flu) = Yas). 


Therefore the result of Section H gives immediately 


2g) L572 (7 OF 


If we have the special case of orthogonal curvilinear coordinates, then of course 





J? = 911922 +--+ IJnn; 
1 


“9 
Gii 


a __ 


so we have 





Vf = a (Lom OF 


Gi Ou; 





For spherical coordinates on R® this gives 
5 1 O a; r siniwoT* 0 (r’sing Of* 
VE = r sin y | = 
r2siny | Or Or + r2 Ow 060 \ r? sin* py 00 
1.@ f sof" 7 1 O of* 1 off 
= i i * 
ro Or Or r? sin yp Op a Op r2 sin? » 06? 
PROBLEM 15-19. For the case of cylindrical coordinates in R*, show that 


vai? oy” | Lory" | oe 
Or Or "72 902° Qz2° 


























L. Parabolic coordinates 





Here’s an example that we have not yet mentioned. It’s a coordinate system for R? based 
on squaring complex numbers. Namely, in complex notation 


1 
rtiy= 5 (u + iv)?. 


24 Chapter 15 


In real notation, 


(: = 3 0%), 


y =U. 


We assume that u > 0, —oo < v < oo. Then we produce all points of R? except the negative 
x-axis (—oo, 0] x {0}. 


PROBLEM 15-21. We call u, v parabolic coordinates because the curves in the x — y 
plane on which wu is constant and also those along which v is constant are parabolas. Prove 
this, and also prove that the origin is the focus of each of these parabolas. 





Our notation from Section E leads us to the frame vectors 


yi(u, v) _ (u,v), 
po(u, v) = (—-v,u). 


These vectors are orthogonal, so we have here orthogonal curvilinear coordinates. 


PROBLEM 15-22. Show that 
1 02 * 02 * 
Vf (0,4) = (S5+54). 





u2 + v2 \ du? ‘3 Ov? 





